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ABSTRACT 

Starting from a Calogero-Sutherland model with hyperbolic interaction confined 
by an external field with Morse potential we construct a Heisenberg spin chain 
with exchange interaction oc 1/ sinh 2 x on a lattice given in terms of the zeroes of 
Laguerre polynomials. Varying the strength of the Morse potential the Haldane- 
Shastry and harmonic spin chains are reproduced. The spectrum of the models 
in this class is found to be that of a classical one-dimensional Ising chain with 
nonuniform nearest neighbour coupling in a nonuniform magnetic field which allows 
to study the thermodynamics in the limit of infinite chains. 
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Studies of Heisenberg spin chains with exchange coupling proportional to certain functions 
decaying like the inverse square of the distance have unveiled surprisingly simple properties 
which allow for a detailed study of these many-body systems. Among of these properties are 
exact wave functions of Jastrow product form and very simple spectra. Most notably, the 
Haldane-Shastry model [p], @] with trigonometric {pal/ sin 2 x) exchange admits an interpreta- 
tion as the generalization of the concept of an "ideal gas" to the case of fractional statistics 

All of these models can be related to Calogero-Sutherland models of particles with an inter- 
nal degree of freedom (spin) using the exchange operator formalism || . Within this approach, 
the spin chains are obtained by "freezing out" the kinematic degrees of freedom in these models 
thereby giving a pure exchange model of particles interacting with inverse square potential on 
a lattice given by the static equilibrium positions of the Calogero model. Within this formal- 
ism it has been possible to construct a family of operators commuting with the Hamiltonian 
and among themselves for the Haldane-Shastry model and a model with rational exchange 
related to the 1/x 2 Calogero model confined by a harmonic potential well || [| |T(J (due to 
its spectrum this model will be referenced as the "harmonic spin chain" in the following). Al- 
though it has not been possible so far to obtain the Hamiltonian from this construction directly 
this is to be considered as the proof of integrability of these models. 

In the limit of infinite chains the above considerations have been applied to the hyperbolic 
chain with exchange oc 1/ sinh 2 x on a translationally invariant lattice. A generalization of this 
case to finite nonuniform lattices which reduces to the Haldane-Shastry and the harmonic spin 
chain in certain limits is presented in this letter. We start with the Hamiltonian of a iV-particle 
Calogero-Sutherland system of particles interacting with an external field with Morse potential 
given by 
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The classical equilibrium positions of the particles (that will later define the lattice on which 
the spins are placed) are given in terms of the set of nonlinear algebraic equations 

N 
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E ?fr t ?3 + r 2 (zi ~ 1) = V = h ■ ■ ■ , N) , (2) 



where we have introduced Zj = e 2xj . A general solution to (|2]) is unknown. However, following 
the observation by Calogero JO] for 1/x 2 case we assume that the {zj} are the roots of some 
polynomial Pn(z) = Ilf=i ( z ~ z j) obeying the second order differential equation [|l^ 

Ap" N + Bp' N + X N p N = , (3) 



where A and B are polynomials in z and A at is some constant. Choosing Xn = 2Nr we find 
that @ reads 

d 2 p N {y) . ( ,-|, -p\ dp N {y) 
V dy2 +(-y + l + T) — + iV>jv(?/) = , y = 2rz (4) 



with T = 2(r — iV) + 1 from which we identify as the Laguerre polynomial L^(2r^). Let 

•(T) 



us note the following properties of L) 
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1. For T > — 1 all roots of L]y are real positive numbers. 

2. L)y has multiple zeroes if and only if T = —n, n — 2, 3, . . . , N. At these values of T one 
can write 

3. As T = — n + e, £ — > the n roots of L N approaching have the following asymptotic 
behaviour: 

Zj ~ const|£:|~ exp ( — — ^ , j — 1, . . . , n . (5) 
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Introducing variables 2j = 1 + t~zQ and taking the limit of large r the rational limit of the 
Calogero model related to the harmonic spin chain is obtained. In this case the lattice points 
(j are the roots of the Hermite polynomial Hn{()- For the special values of T being a negative 
integer > —N the lattice separates into two parts not coupled by the interaction term in the 
Hamiltonian. Upon rescaling one of them is given by 77-th roots of unity (^), which results into 
a lattice of equally spaced sites with periodic boundary conditions. The corresponding spin 
chain is the Haldane-Shastry model with interaction oc 1/ sin 2 x. The other part corresponds 
to the hyperbolic model with T replaced by — r. Finally, as T — > —N the model reduces to the 
trigonometric model. 

Following Fowler and Minahan M we now consider the a model where N bosonic particles 
are sitting on different points Zj satisfying (Q) and allow the exchange of particle positions 
as the only dynamical process. Denoting the corresponding Hermitian exchange operator for 
particles i and j by Aiij we choose the Hamiltonian of this system to be 

N ZZk 
'Hex = X! hjk-Mjk , hj k = - — — . (6) 
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We now construct operators 
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hj = -(ifj + ujj)(itj — ujj) — r (zj — 1) + const ujj (7) 
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where 

N z ■ + z k N 

ftj = E j _ * M ik , = E • 

At some value of the constant in (|7p these operators can be factorized: 

hj = CLjOj, af = - (7Tj ± Vj) ± 5t(zj — 1) , 5 2 = 1 . (8) 

Commutators of two operators from this set are a bit more complicated than in the case of 
harmonic chain, 

[hj, h k ] = 2Sr(hjA4jk — Adjkhj) + UjA4j k hj — hjM.j k ujj . 
Using the properties 

MijZi = zjMij , MijZ k = z k Mij for i ^ k ^ j 

M ijk = MikMij = MijMjk = MjkMik 

of the exchange operators M.y one finds that the commutators [hj,H. ex ] = for an arbitrary 
set of solutions to (||) and hence the operators 

TV 

Zm = E K ( 9 ) 

all commute with the Hamiltonian For the Haldane-Shastry model and the harmonic 
spin chain || it can be shown that operators similar to these X m also commute among them- 
selves. We have not been able to prove this in the general case considered here. After tedious 
calculations we were only able to show that [I n ,I m ] = for 1 < n, m < 3. 

As usual the operators Aijk can be related to exchange operators acting not on the particle 
positions but on their internal degree of freedom as they are defined in the space of bosonic, 
i.e. totally symmetric wavef unctions. This leads us to consider the Hamiltonian 
N n a — 1 N 1 

H = ~ E h 3k^~^ = C N - £ h jk Vj k , C N = —N(N - 1) (3r + 2N-1) (10) 

j<k ~ j<k 

with the spin exchange operator Vjk — % ipj • &k ~ !)• The a" are Pauli matrices acting on the 
j site of the lattice. The same procedure gives a trivial exchange operator X\ = —^J2'Pijk 
for the spin analogue of X±. The next operator in this sequence, Xi, already contains four-spin 
exchange terms with nonuniform exchange couplings. 

To construct eigenstates of the Hamiltonian ( llTf ) we start from the ferromagnetic vacuum 
| ) = | tt ••• T) and consider states with given magnetization 

N AI 
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h<—<3M s=1 
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Using the fact that the Zj are roots of the iV-th Laguerre polynomial we find that the eigenstates 
in the one-magnon sector (M = 1) have amplitudes 



j (r+2m) (l) > 

^(j) oc zT N JrTr , ' m = 0,...,J\T-l. (12) 
Their energies E$ = e m are given by the following quadratic dispersion law: 



TYl 

e m = -{T + m) . (13) 

Next we have studied the two-magnon sector (M = 2). The amplitudes ip(ji, jz) can be 
written as polynomials in {^J -1 }. No closed expression for these amplitudes has been found. 
Nevertheless, the eigenproblem can be solved analytically and we were able to find the complete 
set of iV(iV — l)/2 eigenvalues which can be written as 

E^ n = e m + e n (l-5 m , n . 1 ) , < m < n < N - 1 (14) 

with the single magnon energies e m given by (|TB|). 

Finally we considered the M- magnon sector. Within the Ansatz 

,0 _ n\>n( z j x ~ z j l i) 2 F( z h--- z 3M) 

Y\u=l z jvP'N(. z j„) 

where F is some symmetric polynomial in {z}, the corresponding part of spectrum which 
comprises m^n-^2M+i)\ e ig enva hies is found to be additive 

M 
k=l 

with the dispersion (|T3|) where the integers obey the set of restrictions 

< mi < m 2 - 1 < . . . < m M - 1 < N - 2 . (17) 

While our solution of the M-particle sector is not complete we make the following hypothesis 
concerning the spectrum of the class of models (PH|): 

All the eigenvalues of H can be written in compact form 

N-l 

E ni ...n N = E e fc^+i(! - n k) , (18) 
fc=l 

where = \k (Y + k) and {rik} = 0, 1. As a consequence, the operator H = —1/2 Y,j<k hjk&j&k 
is equivalent up to unitary transformation to the Hamiltonian of classical ID Ising chain in a 
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non-uniform magnetic field 



N-l 



fe=l 



N-2 



— CN-l^N + X! Wk+l^k — Cfc+l) — CTfe+lCit+2 £ it+l] 
k=0 

with {crfc} = ±1. 

We have checked this hypothesis numerically to give the correct spectrum for lattices of 
lengths up to N = 12. We have also confirmed it in the limit r —>■ oo by establishing the 
analytic correspondence to the effective Hamiltonian of the harmonic chain (Eq. (20) in ||). 
Finally, it is consistent with the picture of lattice separation at —N < V = —n < — 1 (where 
e n = 0) mentioned above and coincides with the Ising Hamiltonian that follows for the Haldane- 
Shastry model in the appropriate limit of the Hubbard model with 1/r hopping |13[ . 



The simple form of the spectrum ( |18D allows to compute the free energy in the thermody- 
namic limit by using the transfer matrix method. The partition function on a finite lattice can 
be written as 
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where 
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To perform the thermodynamic limit iV — > oo we rescale the magnon energies ([13]) with a factor 
1/N 2 and find the leading term of ( |19|) to be Z ~ exp(f3E ) Yl^S 2 (1 + w k ) with some proper 
E renormalizing the ground state energy to zero. From this expression one obtains for the free 
energy per site 

(21) 



'=4 4/o 1&hl ( 1 + e " <W/2 ) ' 



Using the quasiparticle dispersion e(x) = x(7 + x) (the renormalization of the magnon energies 
leads to 7 = T/N when taking the thermodynamic limit of (|T3|)) this can be written for 
-1 < 7 < as 



/ = _I Qf 7 rfx i n (i + e ^W/2) + J 1 dxln ^ + e -M*)/2 



(22) 
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Note that changing the sign of the quasiparticle dispersion or, equivalently, the exchange in (|T0|) 
from ferromagnetic to antiferromagnetic changes the free energy by a temperature independent 
term (due to the different ground state energy) only. This has been noticed in the case of the 
Haldane-Shastry model 7 = — 1 before || and is a consequence of the possibility to describe 
these models in terms of an effective Ising model (|18|). 

Since the total spin of the effective Ising model is proportional to the z-component of the 
total spin in the original model one can also study the effect of an external magnetic field. For 
T > — 1 the largest eigenvalue of ([H]) in the M-magnon sector is 



from which we conclude that the ground state magnetization of the anMerromagnetic chain in 



a magnetic field h is given by M.[h) = | (v / 7 2 ~+~8^ — 7 J f° r h < h s = | (7 + 1). Beyond h s 
the magnetization is saturated at |. For finite temperatures the transfer matrix method yields 



In this letter we have constructed a new family of solvable Heisenberg spin chains. As in the 
previously known Haldane-Shastry and harmonic models their spectra can be given in terms 
of an effective Ising Hamiltonian. This mapping allows to study the thermodynamics of this 
models in detail. Generalization to other exchange type models such as interacting SU(N) 
spins or electrons subject to a supersymmetric t—J or Hubbard type Hamiltonian are easily 
constructed by considering different representations of the permutation operators in (fLOf). At 
the same time, a number of interesting questions, especially regarding the construction of a 
complete set of conserved quantities, remains open. 

This work has been partially supported by the Deutsche Forschungsgemeinschaft under 
grant no. Fr 737/2-1. One of the authors (V.I.) has received partial support within the frame- 
work of the Heisenberg-Landau program. 




for M < -N 
~ 2 
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